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Abstract 

The formation of the magnetic moment in C- and N-doped MgO is the result 
of a dehcate interplay between Hund's coupling, hybridization and Jahn- Teller dis- 
tortion. The balance depends on a number of environmental variables including 
electron doping. We investigate such a dependence by self-interaction corrected 
density functional theory and we find that the moment formation is robust with re- 
spect to electron doping. In contrast, the local symmetry around the dopant is more 
fragile and different geometries can be stabilized. Crucially the magnetic moment 
is always extremely localized, making any carrier mediated picture of magnetism in 
(P magnets unlikely. 

In the last few years an intensive research effort has been devoted to synthesize diluted 
magnetic semiconductors (DMS) in the hope of finding a material with semiconducting 
properties and ferromagnetism at room temperature (RT). Standard DMS are produced 
by doping conventional semiconductors with transition metals [1], and the magnetic cou- 
pling usually originates from some carrier-induced mechanism. More recently two main 
experimental facts have challenged this "traditional" picture. Firstly magnetism was 
claimed in oxides DMS on both sides of the metal insulator transition [2], indicating 
that carriers alone are not sufficient to explain the magnetism and that intrinsic defects 
play a crucial role |3]. Secondly, defect-rich or intentionally p-doped oxides revealed the 
possible evidence for RT magnetism. This second class of phenomena has been named 

magnetism, indicating that no ions with partially filled d-shells are at the origin of 



the magnetic moment. The examples of c?" magnets are many and include thin films 
of HfOai], TiOsO El IZ], In203[8], C-doped ZnOP [10], and nanoparticles of different 
materials [TT]. Furthermore studies over chemically synthesized powders suggest that 
magnetism can be found in bulk materials and not only in surfaces [T2| [T3] . 

Usually the formation of the magnetic moment is explained in terms of spin-polarized 
holes localized either at cation vacancies molecular orbitals [TH [T5] or at the p-orbitals of 
the doping impurity [151 [HI [HI [H]. However, a clear understanding of the driving mecha- 
nism behind dP ferromagnetism is still unavailable. In order to understand its uniqueness 
we recall that the record Curie temperature (Tc) for GaAs:Mn (the DMS prototype) is 
200 K[T9], obtained for 8 % Mn doping and the largest hole concentration achievable. 
Since cation vacancies in oxides hardly reach concentrations exceeding 1 % and both C 
and N cannot be doped abundantly [10], one can conclude that a similar Tc in (i° magnets 
would require a magnetic interaction around 10 times stronger than that between Mn in 
GaAs:Mn. Considering that our argument neglects any considerations about percolation 
and the fact that even the origin of the magnetic moment is not established with certainty, 
it is fair to say that any claim of ferromagnetism should be considered exceptional. 

In this arena almost all the theoretical predictions are based on density functional 
theory (DFT) using either the local spin density (LSDA) or the generalized gradient 
approximation (GGA). These describe the ground state of d^ magnets as ferromagnetic 
and metallic, and in fact most of the time as half metallic. However strong electron 
correlations may play a fundamental role in the magnetic moment formation, which subtly 
depends on the interplay between covalency, Hund's coupling and polaronic distortion 
around the impurity [15]. Corrections to the LSDA/GGA such as the LDA+t/ or the 
self-interaction correction (SIC) schemes often return an insulating ground state and no 
long range magnetism. For instance in rock-salt oxides (MgO, CaO and SrO) doped with 
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N substituting for O [T^ [T71 [TS] , one finds that the extra hole entirely localizes around 
one of the 2p orbitals as a consequence of the coupling with phononic modes. Thus the 
physics of these materials reminds that of the manganites [20] , with the difference that in 
dP magnets the moment is associated to the 2p atomic shell. 

In analogy with the manganites, one expects that the interplay between charge, spin 
and orbital degrees of freedom may lead to a number of cooperative physical phenomena |2Uj 
In particular the long range magnetic order in dP magnets is intimately related to the 
formation of the moment itself [16], with Stoner and spin-wave excitations probably com- 
peting. Thus both the moment and the magnetic coupling becomes sensitive to environ- 
mental variables such as doping, charge fluctuations and temperature. In this work we 
investigate the effects of one of these variables, namely electron doping, over the magnetic 
moment formation in C- and N-doped MgO. In particular we will answer two fundamental 
questions: 1) does the magnetic moment survive to electron doping? and 2) what is the 
response of the lattice to such an electron doping? 

Our calculations are performed by using a development version of the DFT code 
SIESTa[21|, implementing the atomic SIC scheme (ASIC)[22|. The core electrons are 
treated with norm- conserving TrouUier-Martin pseudopotentials and the valence charge 
density is expanded over a numerical orbital basis set, including double-^ and polarized 
functions [21] . The real space grid has an equivalent cutoff larger than 800 Ry. Calcula- 
tions are performed with supercells of 96 atoms including /c-point sampling over at least 
25 points in the Brillouin zone. Atomic coordinates are relaxed by conjugate gradient 
until the forces are smaller than 0.01 eV/A. 

As already mentioned when N or C replace O (Nq and Co) the MgO local cubic 
symmetry is reduced[l5]. In both cases the bonds with Mg are longer than the Mg-0 one 
as a consequence of the different ionic radii. However, for Nq two of the three p-bonds 



3 



contract and the hole locahzes around the remaining long one (Nq configuration), while 
for Co only one bond contracts and the two holes localize around the remaining two long 
bonds (Co configuration). In both cases the system is insulating with either one (Cq) 
or two (No) fully filled p orbitals. We now consider Co and add a fractional charge An 
to the supercell, with charge neutrality ensured by a compensating positive background. 
One then expects that the electronic structure and the relaxation becomes progressively 
similar to that of No. 

Our results are presented in Tab. [H where we list the Mg-C bond lengths, dug-C: and 
the cell magnetic moment, function of An. Importantly all the structural calcu- 
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Table I: Mg-C bond length, dug-c, and supercell magnetic moment /i, of Co as a function 
of the electron doping An. We indicate as [Co-Rel] the relaxed structure obtained from 
an atomic relaxation initiated at the An=0 Co geometry, and as [No-Rel] that initiated at 
the An=0 No geometry. In the bracket beside the bond length we indicate the degeneracy 
of the particular bond. 

lations turned out to be sensibly dependent on the initial conditions for the relaxation. 
In particular we find that all the relaxations initialized at the An=0 Co configuration 
converged to a local geometry presenting 4 long bonds and 2 short ones, i.e. still pre- 



4 



senting the symmetry of Cq- In contrast those initiahzed at the An=0 Nq configuration 
converged to a local geometry with two long, two short and two intermediate bonds. 

The changes in the electronic structure as a function of electron doping are illus- 
trated in the cartoon of Fig. [1] for the two obtained geometries. In both cases the ad- 
ditional charge remains localized at the dopant site and the magnetic moment varies as 
H = {2 — An)/iB. The main difference between the electronic structures of the two ge- 
ometries originates from their different orbital occupation. For the Co-relaxed structure 
the fractional charge progressively occupies the minority empty doublet associated to the 
four long bonds, spreading evenly among them [see Fig. [11(a)] . In contrast the p-orbitals 
of the No-relaxed structure form a set of closely spaced singlets. Hence the additional 
fractional charge occupies the first of the available empty singlets and localized further 
along the direction of the bond of intermediate length [see Fig. [U^b)]. These differences 
persist up to An=l where the Co-relaxed structure is metallic, while the No-relaxed is 
insulating. 

Unfortunately the ASIC functional, although constructs a valuable potential so that 
atomic relaxation can be carried out, does not provide accurate total energies. Therefore 
we cannot distinguish energetically between the two geometries found. We have performed 
additional ASIC total energy calculations at the un-relaxed Co and No geometries, which 
seems to suggest a crossover between the two at a doping of An=0.5, however also these 
calculations are affected by an intrinsic lack of accuracy and therefore the result must be 
taken with caution. 

Given this uncertainty we decided to take a look at the theory for the local vibronic 
coupling. For manganites the Jahn- Teller distortion is driven by the coupling of doubly 
degenerate rf-shell eg states to the only two normal modes having the same symmetry. 
In contrast the 2p levels of the N in MgO form a set of degenerate tiu levels, which 
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Figure 1: (Color on line) Cartoon representing the level occupation of Co upon electron 
doping. In panel (a) the energy levels are those of the local geometry obtained from a 
relaxation initiated at the Co configuration, while in (b) those obtained from a relax- 
ation initiated at the Nq configuration. For An=0, in both cases a doubly occupied 
spin-polarized doublet (d^) is followed by a doubly occupied unpolarized singlet (s) and 
separated by the Fermi level (red dashed line) from its minority spin component (d^^). 
Upon doping [An > 0), in (a) the minority doublet gets progressively filled, while in (b) 
the doublet further splits into two singlets and the additional fractional charge distributes 
over the lower of the empty minority singlets. In both cases the magnetic moment goes 
as /i = (2 — An)/iB. 

are coupled with normal modes of symmetry Cg and t2g. Importantly, in this case the 
eg modes tend to stabilize a tetragonal distortion, while the t2g a trigonal one resulting 
in a competition between the two[23j. Furthermore, when anaharmonic corrections are 
considered, a structure of C2 symmetry minimizes the energy. It is notable that the 
highest occupied level for one hole over the impurity in this geometry is 1/ V^{px ± 
Py) and could account for the orbital occupation found for the Co-relaxed structure at 
An=l. Unfortunately, our relaxations always ends up with a tetragonal distortion and the 
symmetry around the impurity is never of C2-type. The vibronic theory then predicts that 
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the ground state is the one with the hole locahzed on just one of the three p orbitals. Thus 
we conclude that simple symmetry arguments seem to support the relaxation initiated by 
the No atomic coordinates, i.e. the one converging to a local geometry having two long, 
two short and two intermediate bonds. 

After having studied Co, we briefly take a look at No under doping. This time we 
find that the unpaired singlet gradually fills upon electron doping, thus that the distortion 
gradually reduces to a perfectly cubic symmetry and the moment follows /i = (1 — An)/iB. 
Interestingly we find that the residual hole always localizes over the longer of the Mg-N 
bonds and that no magnetic coupling is found for every An, despite the material remains 
metallic. 

In conclusion, we investigated the effects of electron doping on the magnetic moment 
formation of N- and C-doped MgO. We find that for all doping concentrations the impu- 
rity levels are deep in the MgO gap and the magnetic moment is stable against charge 
fluctuations. However, lattice distortion always promotes the localization of the doping 
hole, thus reducing the chance of long range magnetic coupling between impurities. 

This work is sponsored by Science Foundation of Ireland. Computational resources 
have been provided by the HEA IITAC project managed by TCHPC and by ICHEC. 
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